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Abstract 

This work derives sufficient conditions for the permanence and ergodicity of a 
stochastic predator-prey model with Beddington-DeAngelis functional response. The 
conditions obtained in fact are very close to the necessary conditions. Both non¬ 
degenerate and degenerate diffusions are considered. One of the distinctive features 
of our results is that our results enables characterization of the support of a unique 
invariant probability measure. It proves the convergence in total variation norm of the 
transition probability to the invariant measure. Comparisons to existing literature and 
related matters to other stochastic predator-prey models are also given. 
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1 Introduction 


This paper focuses on stochastic predator-prey models with Beddington-DeAngelsis func¬ 
tional response. In ecology, a functional response is the intake rate of a consumer as a 
function of food density. It is associated with the numerical response that is the repro¬ 
duction rate of a consumer as a function of food density. Holling [6] initiated the study 
of functional response, where he introduced several types of such responses. The so-called 
Holling type II functional response is characterized by a decelerating intake rate following 
from the assumption that the consumer is limited by its capacity to process food. Similar 
to Holling-type functional response with an extra term describing mutual interference by 
predators, Beddington [Tj and DeAngelis et. al. |3] introduced the nowadays well-known 
Beddington-DeAngelis functional response; see also |24J and the references therein. Such a 
model represents most of the qualitative features of the ratio-dependent models but avoids 
the “low densities problem.” 

As the building blocks of the bio- and eco-systems, the basic premise of the predator-prey 
models is that species compete, evolve, and disperse for the purpose of seeking resources to 
sustain their struggle and existence. Denote the two population sizes at time t by x(t ) and 
y{t), respectively. Then a general deterministic model called Kolmogorov’s predator-prey 
model takes the form 

f x(t)=xf(x,y), 

\ y(t) =yg(x,y). 

When f(x, y) = b — py and g(x, y) = cx — d, one gets the so-called Lotka-Volterra model. 

In addition to the study of deterministic models, stochastic predator-prey models have 
received increasing and resurgent attention. Stochastic models can be considered as the 
above systems subject to Brownian motion perturbations. Rudnicki 123 provided a detailed 
analysis for stability in distribution of a stochastic Lotka-Volterra model. Meanwhile, Mao 
et. al. jTH] and Du and Sam [4] studied general stochastic Lotka-Volterra models using 
Lyapunov-type functions and exponential martingale inequalities. Recently, Lotka-Volterra 
models in random environment have also gained much attention 1251 . In addition, there is a 
resurgent interests in treating evolutionary games [5], in which Lotka-Volterra type equations 
are one of the central models. Concerning different functional responses, references [T6] and 
HZ! dealt with the stochastic predator-prey model with Holling functional response of the 
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form 


{ dx(t) = x(t) (ai — bix(t) — Ciy(t) )dt + ax(t)dBi(t), 

+ *%) (i-i) 

dj/(i) = J/(*)( - - b 2 y(t) + x + x ( t ) ) dt + 

where a*, b i: c. a, and f3 are appropriate constants, and Hj(-) are standard Brownian mo¬ 
tions. Ji et. al. [10] studied the predator-prey model with modified Leslic-Gower and 
Holling type If schemes with stochastic perturbation; see also [9] in which stochastic ratio- 
dependent predator-prey models were considered. Moreover, several stochastic models with 
the well-known Beddington-DeAngelsis functional response were also studied in [HI fl5l [23]. 
In ecology models, an important concept is stochastic permanence, which indicates that the 
species will survive forever. Much effort has been devoted to finding conditions needed for 
stochastic permanence. In some of the aforementioned papers, using suitable Lyapunov-type 
functions, some conditions for extinction or permanence were also provided and ergodicity 
was investigated; see pin]. However, as shown later in Section [4] of this paper, their con¬ 
ditions are restrictive and not close to a necessary condition. In other words, there is a 
considerably large set of parameters satisfying neither their conditions for extinction nor for 
permanence. Moreover, their results are not applicable to degenerate cases. Thus, although 
interesting, their work left a sizable gap. One of the main goals of this paper is to close this 
gap. We aim to providing a sufficient and almost necessary condition for permanence (as 
well as ergodicity) for the following model with Beddington-DeAnglesis functional response, 


dx(t) = x(t)(ai — bix{t) — 


ciy(t) 


dy(t) = y(t)(- a 2 - b 2 y{t ) + 


rrii + m 2 x(t) + m 3 y(t) 
c 2 x(t ) 


nii + m 2 x(t ) + m 3 y(t ) 


)dt + ax(t)dBi(t), 


)dt + 0y(t)dB 2 (t), 


( 1 . 2 ) 


where are positive constants for i = 1,2, m 3 > 0, a ^ 0, (3 ^ 0, and Bi(-), B 2 (-) 

are two mutually independent Brownian motions. When m 3 = 0, the functional response is 
said to be of Holling type II. Moreover, in this paper, we also consider the degenerate case 
Bi(-) = ft (-)• 

The rest of the paper is arranged as follows. Section [2] derives a threshold that is used 
to determine extinction and permanence. To establish the desired result, after considering 
the dynamics on the boundary, we obtain a threshold A that enables us to determine the 
asymptotic behavior of the solution. In particular, it is shown that if A < 0, the predator 
will eventually die out. In case A > 0, the solution converges to a stationary distribution 
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in total variation norm. Moreover, ergodicity is established. Section [2] concentrates on 
non-degenerate case, whereas Section [3] treats the degenerate case B i(-) = B 2 {-). In the 
degenerate case, under usual conditions imposed on the Lie algebra generated by the drift 
and the diffusion coefficients, we investigate the controllability of the associated control 
systems and used certain results in [14] to prove analogous results to the nondegenerate 
case, namely, the existence and uniqueness of an invariant probability measure as well as 
the convergence in total variation of the transition probability. Moreover, the support of the 
invariant measure is described. Finally, Section [4] provides further discussion and insight. 
Among other things, it points out that the techniques used in this paper can be applied to 
other stochastic predator-prey models. 


2 Threshold Between Extinction and Permanence 


Let (fl, T, {-Fi}t>o, P) be a complete filtered probability space with the filtration 
satisfying the usual condition, i.e., it is increasing and right continuous while J- o contains 
all P—null sets. Let B\{t) and B 2 {t) be two JVadapted, mutually independent Brownian 
motions. It is well known that for any initial value (x(0),?/(0)) G M+° (the interior of K+), 
there exists a unique global solution to (ll.2j) that remains in almost surely (see [Ej). To 
proceed, we first consider the equation on the boundary, 


d(p(t) = (p(t)(ai — bi<p(t))dt + otip{t)dB\{t). 


( 2 . 1 ) 


By comparison theorem, it is easy to check that x(t) < tp(t) Vt > 0 a.s. provided that 
x(0) = <£>(0) > 0 and y(0) > 0. If oq < ck 2 /2, we can easily verify item (2) of [T| Theorem 3.1, 
p. 447] to show that lim ip(t) = 0 a.s. Hence, lim x(t) = 0 a.s., which implies lim y(t) = 0 

£—>•00 t —>-oo £—>■ oo 

almost surely (a.s.). For this reason, in the sequel, we suppose that oq > a 2 /2 throughout 
the rest of the paper. 

Defining 6(t) = ln^(t), equation (12.ip becomes 

dd{t) = ^ai ——— b\ exp ( 9(t))jdt + adB\{t). (2.2) 


By solving the Fokker-Pla.nck equation, it is shown that the process 6(t) has a unique station¬ 


ary distribution with density given by f*(x) = C exp ( qx-aexp(x )), where q = 


a = 


2bi 


2aq 

CK 2 


1 > 0 , 


OL 


> 0, and C is the normalizing constant. Since 6(t) = In <p(t), it can be easily seen 
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that <p(t) has a unique stationary distribution //_(•) with density (j>*{x) = Cx q ~ l e~ ax , x > 0. 
It turns out that that C = a q /T(q) with T(-) being the Gamma function and that //_(•) is 
the Gamma distribution with parameters q and a. 

By the strong law of large number type result {21, Theorem 3.16, p. 46], we deduce that 


Jim \ [ tp p (s)ds = ^ ^ [ x p+q l e ax dx = ^ := K p < oo a.s Vp > 0. 


t—>oo t 


r (?) 


chT(g) 


(2.3) 


In particular, with p = 1, A'i = — = —-- ^ . This property implies that 

a b i 


lim -ln<p(t) = lim ( - / faq-6i(p(s) N )ds ] + a lim 

t —^oo t t— »oo \t J o \ 2 / J t—too 




= 0. 


Consequently, 


and 


lim sup - lnx(f) < 0, 

t—too t 


lim sup- / x p (s)ds < K p . 

t —> OO t Jo 


Let be the solution to 


(2.4) 

(2.5) 


= ip(t )( ~ Gl + “— b 2 ip(t))dt + ft'ip(i)dB 2 (t ). 


Then y{t) < ij}{t) \/t > 0 a.s. provided p(0) = -0(0) > 0. Hence, with probability 1 

[j]lim sup - In y(t) < 0, 

t—> oo t 


( 2 . 6 ) 


1 In J. Bao and J. Shao, Permanence and Extinction of Regime-Switching Predator-Prey Models, SIAM 
J. Math. Anal., 48(1), 725739, they said that an additional condition needs to be added to obtain (12.61) . 

C2 

namely —02 —-1-> 0 because it is the condition for ip(t) to be an ergodic process. 

2 m 2 

02 

However, m can be obtained without this condition. Indeed, suppose that —a 2 — - -1-< 0. Let 

2 m 2 


Jt(t) be the solution to 


djj(t) = j){t) ( — a\ -\ - —)dt + pij)(t)dB 2 (i). 

m 2 


We have that 


lim 

t—too 


ln'i/Jt) 


= lim 

t —>OO \ 


( ln^i(O) 


-a 2 - 


P 


C2 B 2 (t) 

m 2 


t 


= —a 2 - 


P 2 


C2 

m 2 


< 0 a.s. 


By the comparison theorem, 


lim sup - lny(f) < lim sup - ln^(t) < lim - ln^i(t) < 0 a.s. 

t—>oo t t—too t t^too t 

Thus. (12.61) holds without an additional condition. 


5 










and 


1 P 

limsup - / y p (s)ds < K p for some constant K p > 0. 

t—yoo t Jn 


(2.7) 


Define the threshold 

A .— —do — - "f" 

" 2 


c 2 x . , , (3 2 a q f°° c 2 x q e~ ax 

/ -r- P-{dx) = -a 2 ~— + t / --- 

l 0 m i + m 2 x 2 i (q) J 0 m i + m 2 x 


dx. 


Theorem 2.1. If X < 0, then the predator is eventually extinct, that is, lim y[t) = 0 a.s. 

£—>■00 

Moreover, as t — » oo t/ie distribution of x(t) converges weakly to //_(•) that is the Gamma 

2 ai 2bi 

distribution with parameters q = —-— 1 and a = ——, respectively. 

a 2 a 2 

Proof. Let y(t) be the solution to the equation 


*/(<) = !/(*)( - “2 “ &2t/(<) + mi -‘ lf + 


( 2 . 8 ) 


where is the solution to (12.11) . By comparison theorem, y(t) < y(t ) a.s. given that 
99 ( 0 ) = x(0),|7(0) = ?/(0). In view of the Ito formula and the ergodicity of <p(t), 

1 


lim sup 

t—KX) t 


In y{t) = lim sup f 7 / ( - a 2 - - b 2 y(s) + C p*\ /cA ) ds + P ^ ^ 

t—*oo \t J 0 2 


mi + m 2 <y9(s) 


, 1 P 

< lim - / 

00 t Jo 


a 2 -^ + 

2 mi + m 2 ip(s) 


CMS) )ds + ft lim Ml = X< 0 a.s. 

7 t—1 OO t 

(2.9) 

That is, ?/(t) converges to 0 at an exponential rate almost surely. The remaining part of the 
assertion can be proved by the arguments in [201 Lemma 7]. □ 


Theorem 2.2. If A > 0 , the process (x(t),y(t)) has an invariant probability measure con¬ 


centrated on 


d2,o 


Proof. For any initial value (m(0), 2/(0)) G M+°, we have 


In y{t) = - 


t 


b 2 y(s)ds 


t 


(J 2 c 2 cp(s) 

— a 2 —-—l- 7 7 

2 mi + m 2 ip[s) 


ds 


c 2 <p(s) 


c 2 x(s ) 


tj o Vmi + m 2 ip(s) mi + m 2 x(s) 
1 I 1 / c 2 x(s ) c 2 x(s 


ds 


>1 r 

-tj o 


t ,/ 0 Vmi + m 2 i(s) mi + m 2 x(s) + msy(s) 
c 2 (p(s) 


ds T /3 


B2O) 


p 2 , 

— — — T 

2 mi + m 2 <p(s) 


ds 


if 

t Jo 


(rrq (^( S ) ~~ + + &2 M S ) V s + P^p-- 


mim 2 


( 2 . 10 ) 
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Letting t —> oo, 


lim inf 


and (I2.10p yield that 
1 f* 


t—>oo t J q \TTli 

Similarly, we have 

1 , , 1 rt , .,2 


f—((f(s) - x(s)) + ( C2m3 + b 2 )y(s))ds > A a.s. (2.11) 

Vm ^ v mim 2 / 


t 


In x(t)=- (ai - y 


t 

>~ f 




ci 2 /(s) 


mi + m 2 x(s) + rri 3 y(s) 


jds + ct 


Si(f) 


— f>i</?(s))(is + J J ^(^(s) - x(s)) - + a 


fliM 


It follows from (12. 3 1) . (12.4p . and (I2.12p that 

( — bi((f(s) — x(s)) H-—> 0 a.s. 

V m i / 


1 

lim inf - 

£—>■00 t 


( 2 . 12 ) 


(2.13) 


c 2 


Dividing both sides of (12. lip and (12.131) by — and &i, respectively, and adding them side 

mi 

by side, we have 


i r 

lim inf - / y(s)ds > 

t^oo t J o 


&im 2 m 2 A 


cic 2 m 2 + feic 2 mim 3 + fei6 2 mfm 2 
For 0 < h <m < H < oo, Holder’s inequality yields that 

1 /■* /i d \ I /f /■* 


= : m > 0 a.s. 


(2.14) 


^ Jo 

which implies that 


^-{y(.s)>h}y(^s^)ds ^ I l{y(s)>h}ds 


t 


t 


y 2 (s)ds , 


lim inf - 

t— 1-00 t 


{i/(«)>ft} ds ^ flim inf - / l{ y ( s )>sH/(s)ds) f lim sup - / y 2 (s)ds 
\ t->0O t Jq J V 1-J.OO t J Q 


-1 


> flim inf- / y(s)ds — h] flim sup- f y 2 (s)ds] > 

t-+ oo f ./ n ) \ t . ^ t In / “ 


(m — ny 

k. 


(2.15) 


a.s. 


In addition, (12.51) and (12.7p imply that 

1 rt ^ i rt 

lim sup - / l{ y ( a )>H}ds <—limsup- / y(s)ds <-£, a.s., 

t—yoo t J q -tl t—>00 t J q -Tj 

1 f* , 1 , 1 /■* . , Ad 

limsup - / l{x(s)>/r}as <—limsup - / x(s)ds < — a.s. 


(2.16) 


£—>•00 


£->-oo 


H 


It follows from (12.15P and (12.16P that for h < —, H > + Ki)K 2 


2 ’ m 2 


lim inf - 

t—400 t 


L{(a:(s),y(s))gn} ( ^'S d 


(m — A ) 2 Ad + Ad m : 


Ad 


> 


8Ad 


a.s., 


(2.17) 
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H 




























(2.18) 


where A = {(x, y) : 0 < x < H, h < y < H}. By virtue of Fatou’s Lemma, we have 
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m 


lirninf- / P(s, (x, y), A)ds > V(x, y) G 
t J o 8 K 2 




where P(t, (x,y), ■) is the transition probability of (x(t),y(t)). By the invariance of Ml — 
{x > 0,y > 0} under equation (11.21) . we can consider the Markov process (x(t),y(t)) on the 
state space A4. It is easy to show that ( x(t),y(t)) has the Feller property. Thus, inequality 
(12.18ft implies that there is an invariant probability measure /i* on Ai] see [I9j. Since 
y(t) —* 0 provided that x(0) = 0, lim^oo P(t, (0, y), K) = 0 for all compact set K C M.. 
Thus, we must have y*({x = 0 ,y > 0}) = 0 (equivalently /i*(M+°) = 1). Furthermore, by 


the invariance of M^°, y* is an invariant probability measure of (x(t),y(t)) on 


d2,o 


□ 


Since Bi(-) and B 2 (-) are independent, the diffusion is non-degenerate. It is well known 
that the existence of an invariant probability measure is equivalent to positive recurrence. 
Hence, the invariant probability is unique and the strong law of large numbers holds; see 
Theorems 3.1, 3.3]. We have the following result. 


Theorem 2.3. If A > 0, (jl.2[) has a unique invariant probability measure /x* with support 


o2,o 


. Moreover, 

a) For any y*-integrable f(x,y ) : M+° —* M, we have 


lim - / f(x(s),y(s))ds= / f(x,y)fj,*(dx,dy)a.s.V(x(0),y(0))e 
t^oo t Jq J 


p2,o 


(b) lim ||P(£, (x,y), •) — /!*(•) || = 0 \/(x,y) G M^° where || • || is the total variation norm. 


t—>o o 


Proof. Assertion (a) was proved in |T3j, Theorem 3.3]; we refer to [12, Proposition 5.1] or (2j 
for the proof of assertion (b). □ 


As a direct corollary of Theorem 12.31 if A > 0, system (1 1.2 j) is stochastically permanent 
in the sense that for any e > 0, there is some 5 G (0,1) such that lim inf P(t, x, y, [A, 5" 1 ] 2 ) > 


£—>•00 


1 — e. Moreover, it follows from (12. 5 j) and (12.71) that we have the following limits. 


lim - / x p (s)ds = I x p y*(dx,dy) a.s. V(x(0), y(0)) G M+°, p > 0, 


t—>oo t 


lim - / y p (s)ds = I y p y*(dx,dy) a.s. V(x(0),2/(0)) G M+°, p> 0. 
t^oo t L I 















3 Degenerate Case 


Suppose that B i(-) = B 2 (-) = kP(-). We consider the system of equations 

ciy(t) 


dx(t ) = x(t)(ai — b\x(t) 


dy(t) = y(t)(- a 2 - b 2 y(t) + 


mi + m 2 x(t) + m 3 y(t ) 
c 2 x(t) 


)dt + ax{t)dW (£), 


-)dt + /3y(t)dW(t). 


(3.1) 


+ m 2 x(t ) + m 3 ?/(f)' 

Owing to the symmetry of the Brownian motion, we can suppose a > 0. Since estimates 
in the previous section still hold for this case, we have lim y(t) = 0 when A < 0 while x(t) 

£—>■00 

converges weakly to the stationary distribution /i_ of <p(t). In what follows, we suppose A > 0 
for which the process has an invariant probability measure y* on M^°. Putting £(i) = lnx(t) 
and rj(t) = In y(t), equation (13.111 becomes 


cv 

d£(t) = («1 - IT - h ^ m -- m , „ (t) 

2 mi + m 2 e^> + m 3 e ??w 


Cie 


v(t) 


B 2 “ “ r „ P £(*) 

Mt) = ( - a 2 - — - b 2 e +-—-77 w —--7 T) 

v 2 mi + m 2 e^w + m^w 


)dt + adW (£), 
)dt + pdW(t). 


(3.2) 


Denote by (£ u ’ v (t),r] u ’ v (t)) the solution with initial value (u,v) to (I3.2[) and let P(t, (u,v), •) 
be its transition probability. Put 


A(u,v ) = 


/ 


V 


a i 


cr 

~2 


feie 1 ' 


Cie 


/3 2 


mi + m 2 e“ + m^e v 
c 2 e u 


—a 2 - b 2 e v H- 

2 mi + m 2 e u + mae^ 


\ 


/ 


and B(u, v ) = 


a 

0 


To proceed, we first recall the notion of Lie bracket. If X(x) = (X!,X 2 ) T and Y(x) = 
(Pi, Y 2 ) T are vector helds on M 2 then the Lie bracket [A", Y] is a vector field given by 


dY, 


dXi 


1 X > Y U*) = [ X ^)- Y , dxi 


dY. 


OX, 


(x) + ( X 2 — —(x) - Y 2 — —(x) , i = 1,2. 


' dx' 


dx 2 


We impose the following condition. 


Assumption 3.1. The Lie algebra L(u,v) generated by A(u,v), B(u,v) satisfies dimL(u,v) = 
2 at every (u, v ) el 2 . In other words, the set of vectors A, B, [A, B], [A, [A, B]\, [B , [A, B]\,... 
spans M 2 . 


This assumption appears to be satisfied for most practical situations. It seems to be 
satisfied for any oq, bi, Ci, mi, m 2 , m 3 , a > 0, i — 1, 2, ft ^ 0 and cq — a 2 /2 > 0, although 
verifying this assumption for our model in general involves cumbersome calculations. For 
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specific parameters, the assumption can be verified by direct calculations. Note that the set 
of (u, v ) at which vectors A, B , [A, B], [A, [A, B]\, [B , [A, B]\, ... do not span M 2 is roots of 
a system of equations det(A,J3) = 0,det(A, [A, B]) = 0,... each of which is a polynomial 
equation of unknowns e u ,e v . Thus, we can show that there is no (u,v) satisfying the above 
system of equations after taking into account a sufficient number of these equations. 

To describe the support of the invariant measure /i* and to prove the ergodicity of 
we need to investigate the following control system 


iiAt) = adit) + ai — —- b\e w AB - 

2 mi + rri2e u A t ) + m$e v A t ) 

B 2 c 

= PW) - fl 2 - -y - b 2 e v AB +-—- (} --tw, 

2 m i + m 2 e ^ ’ + m 3 e^w 


Cie 


y>(d 


(3,3) 


where (j) is taken from the set of piecewise continuous real valued functions defined on R + . 
Let (u^tyUjv), v^tyUjV)) be the solution to Equation (I3.3jl with control 0 and initial value 
(u,v). Denote by Of(u,v) the reachable set from (u,v), that is the set of (u',v') G M 2 such 
that there exists a t > 0 and a control </>(•) satisfying = u',u^(f,u,u) = v'. It 

should be noted that Assumption 13. 1 1 guarantees the accessibility of ( 13. 3 j) . i.e., Of(u,v) has 
non-empty interior for every (u, v ) G M 2 (see in])- We first recall some concepts introduced 
in na. Let U be a subset of M 2 satisfying the property that for any w±,w 2 G U, we have 
w 2 G Of(wi). Then there is a unique maximal set V D U such that this property still holds 
for V. Such V is called a control set. A control set C is said to be invariant if Of(w) C C 
for all w G C. 

Putting — v,j) — -u^ we have an equivalent system 


where 


and 


f = a<i>(t) + g{u^,(t), z^t)), 
\ z^t) = hiu^it), z^(t)), 


a 


g(u, z) = ai —-— b\e 


C\e e<* 


m i + m 2 e u + m-ie z e<* r 


(3.4) 


h(u , z) — — (a 2 + — + — (ai ——)) — b 2 e z e aU + —b\e u + 
\ 2 ol 2 / OL 


P t 


c 2 e u + — c\e z+ ^ u 
a 


mi + m 2 e u + m 3 e z+ a u 

Denote by (u, z) the set of (u', z r ) G M 2 such that there is a t > 0 and a control </>(•) such 
that u, z) = u', Z(j,{t , u, v ) = z'. 
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Claim 1. For any Mo, Mi,zo G M and e > 0, there exists a control <fi and some T > 0 such 
that U(p(T , Uq, Zo) = «i, ^(T, wo^o) - z 0 | < £■ 

For the proof, suppose that u 0 < U\ and let p\ = sup{|(/(M, z)\, \h(u, z)\ : u 0 < u < 
Ui, |z — z 0 | < e}. We choose <f(t) = p 2 with (ap 2 pf l — l) e > u± — u 0 . It is easy to check that 
with this control, there is a T G [0, spf 1 ] such that u^T, u 0 , z 0 ) = Mi, \ z^(T, u 0 , z 0 ) — z 0 1 < e. 
If u 0 > Mi, we can construct </>(£) similarly. 

Claim 2. For any Zo > Zi , there is a Mo G i, a control (j), and some T > 0 such that 
z^T, Mo, ^o) = zi and that u^(t, uq, zq) = uq VO < t < T. 

Indeed, if /3 > 0 and — no is sufficiently large, there is a p 3 > 0 such that h(uo, z) < 
—p 3 Vzi < z < zo. This property, combining with (13.4j) . implies the existence of a control (j) 
and a T > 0 satisfying the desired claim. In case f3 < 0, choosing uq to be sufficiently large, 
we have the same result. 


Claim 3. If 0 < (3 < a, for any z 0 < Zi, if u 0 is sufficiently large, inf z e[z 0 , Zl ]h{uo, z) > 
0, which implies that there is a control <f and a T > 0 satisfying z^(T, u 0 , z 0 ) = z i an d 
M 0 (f, m 0 , z 0 ) = m 0 VO < t < T. 

Lemma 3.1. Suppose /3 < 0 or /3 > a. Let c* : = sup < z : sup{h(u, z)} > 0 Vz < z. >. Then 

1 uGR > 

c* > —oo, (c* may be oo) and for any (m, z) G K 2 , Of{u,z) D {(u',z') : z' < c*}. 


Proof. Note that 

In view of Jensen’s inequality, 


A — —cto —-h 

2 


c 2 x 


! 0 mi + m 2 x 


P-(dx) > 0. 


/ 0 mi+m 2 x 


C2X »-(dx)< c >P*»-(dx) 


c 2 (ai -y )^ 1 

o? 

mi +m 2 (ai - y )^ 1 


mi + m 2 J 0 °° xp-(dx) 


_ Q/ 

If e u = (ai — -)bf \we have 


c i(°i - yK 1 


h(u, z) 
Since 


/ \ , _1 t- 

mi + m 2 [a 1 - —) b 1 + m 3 ebm 


ci(ai - %~)b 1 1 


q 2 0 n zJ r — U 

(a 2 + y) + b 2 e z e^ +-^--y • 

z mi + m 2 e u + m 3 e z+ «“ 


- (a 2 + ^r) > 0 , 


mi +m 2 (ai - y )^ 1 


2 
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h(u, z) > 0 when e z is sufficiently small. Now we move to the second assertion. Note 
that it follows directly from the continuous dependence of solutions on initial values that 
if Of(w 2 ) C Of (w\) provided w 2 G Of{w 1 ) (u)i,u} 2 G K 2 ). For (u,z) G R 2 , define = 
sup {zi : 3ui such that (ui, Z\) G Of(u : z)}. For any (u±, Zi) G R 2 , it is easy to derive from 
Claims 1 and 2 that Of(ui, Z\) D {(«', z') : z' < z{\. Hence Of (u, z) D {(u 1 , £ 1 ) : z 1 < 3 U , 2 }. 
If 3«,2 < c*, there is some «Gl such that h(u,$ U:Z ) > 0. Since h(-) is continuous, there is an 
z > 3m ,2 such that inf{h(w, z) : z G \i u ,zi S]} > 0. As a result, there is a control 0 and a T > 0 
such that ^(T,w, 3 m > 2 ) = F and u^u, 3 ^, 2 ) = G [0,T], That is, (u, z) G C>^(w, 3 m, 2 ) C 
which contradicts the definition of $ UtZ . The proof is complete. □ 

Proposition 3.1. The control system f!3.3p has only one invariant control set C. If 0 < 

/3 < a, C = M 2 . If /3 < 0 or (3 > a, C = Uu,v) : v — —u < c*}. 

a 

Proof. If 0 < 0 < a, it follows from Claims 1, 2, and 3 that for any (ui, 2 i), (u 2 ,z 2 ) G M 2 , 

(u 2 ,z 2 ) G Of(u\,Zi). Hence, for any (wi, t>i), (u 2 , v 2 ) G M 2 , we have (u 2 ,v 2 ) G Of(u±,Vi). 

This implies that M 2 is an unique invariant control set. Now, consider the case (3 < 0 

or f3 > a for which the conclusion of this proposition is a direct corollary of Lemma 13.11 

if c* = 00. If c* < 00, it is seen from the definition of c* that h(u,c*) < OVu G M. 

Consequently, for all control 0, we have z^it.u.z) < c*Vt > 0 provided that z < c*. In 

other words, Of(u,z) C {( u',z r ) : z! < c*}. This claim combined with Lemma [3.11 implies 

that Of(u,z) = {( u',z') : z' < c*} for all u G M, z < c*. As a result, {( u',z') : zJ < c*j 

is a invariant control set for (|3.4ji . The uniqueness of this invariant control set is obtained. 

in the property that {( u',z') : z' < c*} C Of(u,z) for every ( u,z ) G M 2 . Equivalently, 

C := {(u,v) \ v — —u < c*} is a unique invariant control set for (I3.3|b □ 

a 

Note that if A > 0, there is an invariant probability measure n* of (13.21) that is associated 
with p* of (13. lj) . Since there is only one invariant control set C , it follows from Assumption 
13. 1 1 that n* is the unique invariant probability measure with support C. Moreover, for all 
(u, u)gC and a 7r*-integrable function / we have 

P{ ( limi /'' f(C’ v (s),v u ’’(s))ds = 

These results are proved in [T3j. Moreover, it follows from [ 121 Proposition 5.1] 

lim ||P(t, (u, v), ■) — 7 r*(•)|| —>■ 0 V(u, v ) G C , (3.6) 

t—>00 

where || ■ || is the total variation norm, if we can verify the following Hormander condition. 


f(u\v')7r*(du',dv') > = 1. 


(3.5) 
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Assumption 3.2. The ideal Lq in L generated by B satisfies dimLo(u,v ) = 2 at every 
(u,v) G C. In other words, the set of vectors B , [A, B], [B , [A, B]\, [B , [B , A, B]\,... spans 
R 2 . 


We aim to prove that (13.5|) (under Assumption 13.IK and (13.61) (under Assumption 13.21) 
hold for all (u , v ) G M 2 . We need only consider the case /3 < 0 or (3 > a since C = R 2 in case 

0 < (3 < a. 

Proposition 3.2. Suppose that /3 > a, X > 0. Then, for each initial value (u,v) G R 2 , we 
have Tqo almost surely with rfA = inf {f > 0 : (£ u,v {t),ri u ’ v (t)) G C°}. 


The proof of this proposition is divided into several lemmas. We consider only the case 
c* < oo since the assertion is trivial if c* = oo. Let us first explain the idea of the proof. 
Denote d\ = In H , d 2 = In h, where h, H are defined as in the proof of Theorem 12.21 Since the 
process is recurrent relative to A := {(w, v) : u < d\, d 2 < v < d\}, in order to show rfff < oo, 
we need to estimate (uniformly) the probability of entering C° from A. The difficulty is that 
A is not compact. Therefore, we divide A into A\ = {(w,u) : u < d§, c ?2 < v < d\} and 
A 2 = A\ Ai, where — d§ is sufficiently large. Noting that A 2 is compact and using the support 
theorem and the Feller property, we can obtain a positive lower bound for the probability of 
entering C from A 2 . To obtain similar result for Ai, we will analyze the property of the drift 
when — u is sufficiently large and then estimate using the exponential martingale inequality. 

Fix 0 < 5 < min{ai — a 2 /2,a 2 + /3 2 /2}. Thus, there is a 4 < d 2 such that for all 
u < a/3 _1 (d3 — c*),v < ds, we have 

.,2 / 3 2 


a 


d\ 


T - i, ' e ” 


cie 

mi + m 2 e u + m 3 e l 


> S and — a 2 


- b 2 e v + - — - 

2 m\ + m 2 e u + m 3 e v 


< - 8 . 

(3.7) 
) < 1 . 


O' §£ 

Let d/i < min{ — (d 3 — c*),d 3 } — I where I > 0 be chosen such that 2exp(— -—— 

p (a + p) z 

Ot 

Construct open sets D = {(u,v) GR 2 : u < —(d 3 — c*),v < d 3 } and E — {(u,v) G R 2 ,u,v < 
d 4 }. Then put E 3 = E° D C°, E 2 = E \ E 3 . 


Lemma 3.2. Suppose that (3 > a. There is a p > 0 such that 

P {C' v (°d V ) = -jf( d 3~ c*), r) u,v (°D V ) < d z} > Pi V (u, v ) G E, 
where ofifi is the first time (f u,v (t), rj u ’ v (t)) exits D. 
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2 ,a 


Proof. Define T uv — — (—(d 3 — c*) — u + £). By the well-known exponential martingale 

0 x p ' 


inequality, we have P(fi 4 ) > p\ := 1 — 2exp(- 


81 


ip. + P) 


-), where 


« + i,y 


For w 6 Hi and u,v < d4 , it follows from the property of fl 4 and (13.21) that 
A T„) >u + 5(a u p A f„,„) - a > y A T„) - 


2 (cr + P) 


ol T P 


5 


(3.8) 


and that 


>u-e + -(a% v AT UiV ), 


r] u,v {aff A T u>v ) <d 4 - d(ap v A T u ,„) + 


P5 


2 (a + /?) 


A T ^) + * < d 3 . (3.9) 


^ 5 ^ ex 

If ap v > T UjV i it follows from (13. 8 p that fi a,v {T u ,v) > u — £ + -T U)V > — (d 3 — c*) which is a 

^ _ P 

contradiction. Hence ofi v < T UjV for all u G Oi. Furthermore, (13.9ft implies that for u 6 Hi, 


a 


r ) u ' v (crp v ) = V u,v ( a L> V A T u ,v) < d 3 and consequently f u ’ v (afi v ) = — (d 3 — c*). As a result, 

P{C’>^) = - 5 (d 3 - c*),r] u,v (op v ) < d 3 } > P(O0 > p^(u,v) E E. 

The lemma is proved. □ 

Lemma 3.3. Suppose that P > a. There are d 3 € M, p 2 > 0 and T > 0 such that 

^{te’ v <T} > p 2 Vn < d 5 , d 2 < v < di, 
where is the first time (£, u,v (t), rj u,v (t)) enters E. 

Proof. It is readily seen that there are G\ < d 4 , G\ > 0 and > 0 such that 


r a , u C l e 

sup \a\ - b\e - 

2 mi + m 2 e u + m 3 e v 




and that 


P 2 , v . c 2 e u 


sup {— a 2 - b 2 e v H- 

u<cti,i>sr 2 mi+m 2 e u + m 3 e l 


}<Gi, 


}<-Si 


Fix S 2 > 0. Define T = 2— - 4 2 and d 5 = cr 4 — d 2 — (G\ 4 —-)T and the stopping time 

01 2 


( u ’ v = inf{t > 0 : C ,v (t) > cri or r/ u ’ v (t) < d 4 }. 
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By the exponential martingale inequality, we have P{fi 2 } > p 2 := 
where 


1 - 


exp( 


SiS 2 

(a + ft ) 2 


Pl 2 



Si 


2 (cr + ft) 


i] < 


S 2 \ 

ol T ft / 


> 0 


For uJ6fi 2 and u < d 5 , d 2 < v < <4, it follows from the property of ff 2 and 03.21) that 


C’ V (C’ V A T) <u + Gi(( u,v AT) + n/ aS : L (( u ’ v A T) + " 2 


2 (a + ft) 


ol + ft 


Si. 


<g?5 + S 2 + (Gi + —)T — a i 7 


(3.10) 


and that 

!)”,’(C“'” A T) A T) + ^ (C" A T) + 

2(a + ft) a + ft (3 11) 

<di + 5 2 - 5 ^{C’ v TT). 

If £ u,,u > T, we deduce from 03. lip that rj u ’ v (T) < d 1 + 5 2 — -2-(T) = d 4 , which contradicts 
the definition of ((“’A Hence for u j G fi 2 , we have (f u,v < T. Moreover, Q3.10p implies that 
< ai j n v j ew 0 f fi ie definition of ( u ’ v , we have rj u,v ((f u,v ) = d 4 in fi 2 , consequently 
r^ 1 '’ < T in fi 2 . As a result, for any u < d 5 ,d 2 < v < <4, P {rff v < T} > P(fi 2 ) > p 2 . □ 


Lemma 3.4. Suppose that ft > a, A > 0. For any (u,v) G M 2 , the process (f u,v (t), r] u,v {t)) 
is recurrent relative to E, that is, there is a sequence of random variables (t n (u;)} such that 
t n {oj) | oo as n —y oo and that (£ u,v (t n ), r) u,v (t n )) G E\/n G N for almost all co. 

Proof. Since Ei C Of(u,v) V (u,v) G M 2 , it follows from the support theorem (see j?, The¬ 
orem 8.1, page 518] or [22]) for diffusion processes, that there is a T U)V > 0 such that 
P {(f u ’ v (T Ut v),i] u ’ v (T UtV )) G -Ed} > 2 p u ’ v > 0. Since the process (f(t),r]{t)) is Feller and E x is 
an open set, there is a neighborhood 14 ,v of (u, v ) such that for P{ (£ u ' ,v '(T UiV ), rj u ' ,v ' (T UiV )} G 
Ei} > p UtV W(u',v') G V U)V . Let d 5 be as in Lemma [3.31 we consider the compact set 
K = {(u,v) : d 5 < u < di, d 2 < v < d{\. By the Heine-Borel theorem, there is a hnite 
number of V UijVi ,i = 1 such that K c Uf =l V Ui ^. Letting T K = ma x{T u . jVi ,i = l,n}, 

p K = min{p u . jVi , i = l,n}. we claim that for any (u,v) G K, P{t£*’ < T K } > P{r](A < 
Tk} 4 Pk > 0. Combining this result with the conclusion of Lemma 13.31 we derive that 
there are T > 0, p > 0 such that 


P( rff v < T) > p\/ (u, v) G A := {(n, v) : u < <4, d 2 < v < d\}. (3-12) 
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Since (12.17ft is equivalent to 

1 0 a.s., V (u, v) G M , 

the process (£, u ’ v (t), r) u,v (t)) is recurrent relative to A. Using this property, the strong Markov 
property and (13.12ft . we can conclude the recurrence relative to E of (£, u ’ v (t), r} u,v (t)). □ 

Proof of Proposition \3.A Since (f u,v (t), rj u,v (t)) is recurrent relative to A and E, we can 
define the following sequences of stopping times 



= inf {t > 0 : e ,v (t),V u,v (t)) e E}, 
v n = inf {t > q n : £ u,v (t), r) u ’ v (t)) G A}, 
Wi = inf> v n : £ u ’ v (t), rj u>v (t)) G E}, 


which are finite almost surely. 

We also define t n = inf {t > q n : £ u,v (t), rj u,v (t)) ^ D}. Since ECDC A c , it is easy to see 


a 


that o n < i n < v n . Consider a sequence of events O n := {£, u ’ v (t n ) = —(d 3 — c*) , rj u,v (t n ) < d 3 }. 
If we are in the time q n , O n is the future information while we have already known whether 
O n _i has happened. Moreover, it follows from Lemma [3721 that = u',r] u ’ v (q n ) = 

v') < 1 — piV(u',v') G E. Hence, using the strong Markovian property of (£“’”(£), rj u,v (t)), 
we can prove that 

n 

P| O^j < (1 — pi) n —> 0 as n —> oo. 

k =1 


This means that almost surely, O n must occur for some n = n(oj). Whenever O n occurs, we 
have (€ u,v (t n ), V u,v ( L n)) C C°. The proof is complete. □ 


For the case (3 < 0, we have a similar result. 

Proposition 3.3. Suppose (3 < 0, A > 0. Then, for each initial data (u,v) G M 2 , rfA < oo 
almost surely. 

Proof. We only consider the case c* < oo for which C = {(u,v) : v < c* — ru } with 
r = — ^ > 0. Let A be as in the proof of Lemma [3.41 Divide A into A\ and A 2 defined by 
Ai = ACC° and A 2 = A\Ai. It is easy to see that A 2 is compact. Using the same arguments 
as in the proof of Lemma 13.41 we can find T^ > 0 such that inf^, v ^ £ x 2 T c° V < ^a 2 ) > 0- 
Since A\ C C°, we have 

inf ^P( Tco V < T A2 ) = inf ^ < T M ) > 0. 

(u',v')gA (u',v')gA2 
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Moreover, since (£, u ’ v (t), r] u,v (t)) is recurrent relative to A, we can use the strong Markov 
property to obtain the desired conclusion. □ 

We complete this section by presenting the following theorem. 

Theorem 3.4. Suppose a, (3 ^ 0, \ > 0, and Assumption \3.1\ holds. Then, (13. 2 p has a 
unique invariant probability measure n* satisfying that for any n* -integrable function f, 

p{ lim - ff(^ u,v {s),rj u,v (s))ds = ( f(u',v')M(du', dv')\ = 1 V(w,u)gR 2 . (3.13) 

l t J 0 V ,/ R 2 ' J 

Moreover, if Assumption [PI is satisfied, the transition probability P(t, (u,v), •) converges to 
7r*(-) in total variation as t —* oo. 

Proof. The assertions can be proved using (13.51) . (13. 6p . Propositions 13.21 and 13.31 □ 


4 Discussion 


We compare our results with some of the recent results in the literature. In j8j Theorem 

c 2 fd 2 

4.1], under the conditions —— < a 2 H-and ai > a 2 / 2, it was proved that the predator 

m 2 2 

will eventually die out while the distribution of x(t) converges weakly to the stationary 
distribution of u(t). In contrast, using Theorem 12.11 of this paper, we obtain the same 
conclusion provided that cq > a 2 /2 and A < 0. Note that A < 0 is equivalent to 

/»r*3 

c 2 x 


A := 


/ 0 mi + m 2 x 


p-(dx) < a 2 + p /2. 


c 2 


It is easy to verify that A < —, which indicates that our result on extinction of predator is 

m 2 

sharper. Furthermore, a suitable Lyapiinov function was used in [8] to obtain the ergodicity 
of system (jl.2[) for the non-degenerate case as follows (see [HI Theorem 3.1]). 


Theorem 4.1. Assume (c 2 — a 2 m 2 )ai/bi > a 2 mi,bi > aim 2 /{m\ +m 2 x*) and a > 0, /3 > 0 
such that S < mm{c 2 (bi—m 2 (ai — bix*)/mi)(mi + msy*)(x*) 2 ,b 2 ci(mi+m 2 x*)(y*) 2 }, where 
5 = c 2 x*a 2 /2 + ciy*/3 2 /2 and (x*,y*) is the equilibrium of the deterministic system 


{ x(t) = x{t){a i — b\x{t )) 
y(t) = {-a 2 - b 2 y(t ) + 


ciy(t) 


mi + m 2 x(t) + m 3 y(t ) 

__ ) dt 

mi + m 2 x{t ) + m^y{t) 


)dt, 


(4.1) 


Then there is a stationary distribution 7r(-) for system (1.2) and it has ergodic property. 
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To show that their assumption is more restrictive than our assumption of ergodicity, let 
G be the space of the positive parameters (a*, 6 *, c*, rrij, a, (3), i = 1 , 2 , j = 1 , 2 , 3, ai > a 2 / 2 , 
and 

G + = {(ai, bi, Ci, rrij, a,j3)\\> 0}, G~ = {(ai, b t , c h rrij, a, (3) : A < 0}. 


It is easy to check that A is a continuous function of parameters. Hence G + and G~ are 
open. Moreover, the closure cl(G”) = {A < 0} = (G + ) c , which is a necessary condition for 
the extinction of the predator. Let J be the set of parameters satisfying the assumption of 
Theorem 14.11 we must have G~ U J = 0. Since J is open, cl(G - ) U J = 0 or equivalently 
J C G+. 


We will show that J is a proper subset of G + . Choose a±, b\, c\, a 2 , c 2 , rrii, i = 1,3 ,a,/3 

such that A > 0. This choice can be done by taking a\ sufficiently large. Now fix these 

parameters. Since A does not depend on fc> 2 , we claim the ergodicity holds for all &2 > 0 . 

It can be proved that there exists M > 0 independent of b 2 such that x*,y* < M, where 

(x*,y*) is the positive equilibrium of (14.1 ft (if it exists). Thus, for sufficiently small b 2 such 

that S > b 2 c\(mi + m 2 x*)(y*) 2 , the assumption of Theorem 14.11 does not hold while A > 0. 

Next we look at the case mi = 1, m 2 = 1, m 3 = 0 for which the functional response is 

said to be Holling type-II (see (11.11) ). We will make a comparison with the findings in [TT j 

o? 

in which they proved that if ai -— > 0 and c 2 + a 2 —— < 0 , the predator will extinct 

while x(t) converges weakly to the stationary distribution of <j)(t). Moreover, it was shown 
that the system is persistent in time-average if 

2 


cr /3 2 , ai — %■ c 2 + a 2 — 

a 1 — — > 0, a 2 —— > (J, and -— > 


2 Ci b 2 

In the same manner as in the previous part, we can show that our conditions for extinction 
or permanence and ergodicity are weaker than those in H7j. 

We have investigated (11.21) and (13.11) when A ^ 0. Note that the set {A = 0} has Lebesgue 
measure zero in the space of parameters G. Although the set {A = 0} is negligible with 
respect to the Lebesgue measure, it is still interesting to explore the asymptotic behavior 
of the solution in this critical case. The question of asymptotic behavior corresponding to 
A = 0 remains open. To treat this case, new techniques are needed. Moreover, it seems 
that our methods are applicable to stochastic predator-prey models with different types of 
functional responses as well as different diffusion coefficients. Furthermore, our method can 
be applied to stochastic models with Markovian switching. 
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